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lapsing star is found. The semiclassical wave function of that model turn out 
to describe a bound state. For fixed initial radius of the collapsing star, the 
^ ■ corresponding Bohr-Sommerfeld quantization condition implies mass quanti- 

zation. An extension of this model, and some consequences, are considered. 
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The simplest model of a spherically symmetric homogenous dust collaps- 
ing star is that of Oppenheimer and Snyder [1]. In that model one connects 
continuously the Schwarzschild interior to the interior region, which is a slice 
of a Friedmann universe. 

To construct the (mini) superspace of that model one needs the relevant 
Lagrangian and Hamiltonian. Generally one does not use the Lagrangian 
description of a dust distribution, because its density p(x) , is not a dynamical 
variable. Instead one uses the energy momentum tensor, = pU^Uy , 
where is the four-velocity of the matter particles. The Einstein equations 
are G^ u = 8tcT^ . In [2] it was shown how one can define a matter dust 
Lagrangian that will lead to those equations of motion. 

We will use geometrical units: G — c — 1 and signature ( — h +- 1-). 
The interior region of the Oppenheimer-Snyder model is described by the 
Friedmann line element, 



The range of x is < x < Xo > where xo < tt/2 . At x — Xo the interior 
is matched to the exterior Schwarzschild solution. At least classically, the 
energy momentum tensor must be conserved, V = . In the case of 
a Friedmann universe, this implies that p = po/a 3 , where po is a constant 
determine by the initial conditions. In the collapse case one gets po = 3ao/87r 
, where ao is the radius of the Friedmann ball at the beginning of the collapse. 
So the density of the star is 



where r s is the surface radius of the collapsing star. 

To get the matter Lagrangian one defines the "weight one density" p = 
yf^gp [2] (which is the generalization of a delta function). Then the matter 




ds 2 = -N 2 {t)dt 2 + a 2 {t)[d X 2 + sin 2 xdtt 2 2 ] 



(1) 



P 



3a /87va 3 (t) , x < Xo (r < r s ) 
, X > Xo (r > r s ) 



(2) 
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Lagrangian is Lm = — 87r / pg tJ,v U IJ ,U v d 3 x . The four-velocity satisfies (classi- 
cally) U^U^ = — 1 , so from (1) and (2) we get 



32n' 



Xo 



sin 2 xdxNa 3 p 



(3) 



The gravitational Lagrangian may be split into its interior and exterior parts, 



J G 



■i-R 



\o 



sin xdx 



3a 



3Na 



+ 



r>r s 



-gRd 3 x . 



(4) 



The total Lagrangian is then 



L G + L 



M 



An 



XO 



sin 2 xdx 



The matching conditions are [1] 

M 
R 



aa 2 



1 



— N(a — a ) 



+ 



r>r„ 



-gRd 3 x (5) 



-a sin xo 
a sinxo 



(6) 



where i?o is the radius of the surface of the star at the beginning of the 
collapse, and M is the mass of the Schwarzschild solution. As one can see 
from (6), xo is determined by the initial condition: 2M/Rq = sin 2 xo ■ For 
a reasonable cosmological initial condition we have Rq » M so: Xo << 1 • 
In that case the integral over the slice of S 3 is 47rxo 3 /3 , and 



M 



1 



RoXo 



Hence 



4ttXo 



aa 



— N(a — a ) 



+ 



r>r a 



-gRd 3 x . 



(7) 
(8) 



The classical solution of (5) (or (8)) is easily obtained^: For r > 



we have the Schwarzschild solution. For r < r s (x < Xo), it is convenient 

2 Note that for our dust, the only variables in (5) (or (8)) are the metric components 
(a, N) since the the matter field p is determined by a; the two corresponding Euler- 
Lagrange equations are of course equivalent to the Einstein ones with the dust source. 
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to choose the gauge for which N(rj) = a(rj). Then the solution for which 
a(r] = 0) = a and 0.(77 = 0) = , is 

a(rj) = a (l + cosr])/2 . (9) 

Our strategy is this: (5) (or (8)) is our model. This is our starting point. 
Classically it describes the Oppenheimer-Snyder collapse process; we will find 
the semiclassical wave function of this universe. This approach is consistent 
with the analysis of [3], in which a general (inhomogeneous) distribution 
of spherically symetric dust matter is considered. The field describing the 
matter (the dust) is eliminated by gauge conditions and one ends up with 
only gravitational degrees of freedom. 

Usually one starts from a known matter distribution and tries to find 
the metric. Here we go in the opposite direction: we first determine the 
geometry (the metric), and then find the matter distribution. Classically 
both directions are equivalent. Semiclassically we can use our model, because 
we need only the classical action. Higher quantum orders will require a 
modification of this modelf]. We will see that the WKB wave function of 
that universe describes a bound state, so we will get quantization conditions. 

One can use the same formulation to study closed Friedmann minisuper- 
space models filled with dust (and/or radiation) [4]. In that case < x < 27T 
, and there is no "outside" Schwarzschild region. 

The semiclassical wave functionF] of this collapsing star is 

i> WKB {t) = A{t)exp[iS C iass.{0,t)/h] (10) 

where 

Sciass.(0, t) = AtxxI f L{t')dt' + / dt' f V=gR S ch W .d 3 x . (11) 

JO JO Jr>r s 

3 Because for example U^U^ 7^ —1, and we cannot get (3). 

4 We work in the Lorentzian section but one can get exaclty the same results working 
in the Euclidian section. This is true also for the standard minisuperspace models, and 
one may argue that it should be the same in the full superspace. 
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Because Rschw. = , only the first term in the r.h.s of (11) will contribute 
to the semiclassical wave function. We can calculate (11) very easily in the 
gauge N(r)) = a(rf) , for which 



C[7rx a (r] - -sin2rj) 



exp (^Kxlaliv - \sin2n) 



(12) 



One can easily see that (12) is the WKB solution of the one-dimensional 
Schro dinger equation 



i/j(a) = 



where 



V(a) = (8irxl) 2 a(a — a ) , < a < a 
The semiclassical solution of (13), (14) is 

^wKB( a ) = C(p(a)y 1/2 exp[ip(a) /h] 

where 

p(a) = p \J E - V(a')da' = f \j\V{a')\da' 

J an J an 



(13) 



(14) 



(15) 



(16) 



2(2a — a )(a a — a 2 ) 1 ^ 2 — Oq (arcsin(\ — 2a/a ) + — 



Using (9), one can see that (12) and (15) are identical. 

The Schrodinger equation (13), (14) is just the Wheeler-DeWitt equation 
of this simple model [5,6,7]. To see this, we go back to (8) and construct 
the Hamiltonian in the standard way. One should remember that only the 
first term in the r.h.s of (8) contributes to the semiclassical dynamics. The 
Hamiltonian is 

H = dP a -L (17) 
where P a = dL/dd = 8iTXoad/N , namely 



H — N 



167TXoa 



Pa + 47r Xo(« - Go) 



(18) 
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The Wheeler-DeWitt equation is the quantum version of the classical Hamil- 
tonian constraint, OH /ON = , in the coordinate representation \if) >= ip(a) 
and P a — > —ihd/da . Using (18) we get 



167rxo 3a 



aa z 



V»(a) = (19) 



where the potential is (14). So (13) is (as we expected) the Wheeler- DeWitt 
equation of this model. 

We see from (14) that there are two turning points, at a = and a = a 
, between which there is a bound statd^ 

Consequently there will be quantization conditions. In our semiclassical 
limit (in which General Relativity should be a good approximation) the Bohr- 
Sommerfeld quantization condition says that 

J sJ\E -V(a)\da= J yJ\V(a)\da = (n + l/2)irh . (20) 

Let Rq be the initial radius of the collapsing star. Then (20) is a quantization 
condition on the mass of the star. From (7), (14) and (20) we get the mass 
quantization condition 

MH = 2^(^) 3(n + 1/2)2Mp (21) 
where lp and Mp are the Planck length and mass (in the geometrical units 

l p = Mp = h l/2 ). 

Before we consider some of the consequences of (21), it is important to 
verify its extension. Is it model dependent, and in what way? 



5 The wave function of this bound state is of course real 
which corresponds to the Hartle-Hawking wave function [7] . 
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In the Appendix we consider in details the addition of a conformally 
invariant scalar field as a small perturbation. The resulting Wheeler-DeWitt 
equation is (see (41), (42)) 

d 2 



-h^ + V m (a) 



where 



3^ 



V m (a) = (8vr^) 



da 2 



a la — do — 



i/j(a) = 



a 



(m + l/2)ti 



(22) 



(23) 



where a is a real number, and m is an integer (it is the quantum number 
of the scalar field inside the star). Comparing (14) with (23), we see that 
small perturbation requires that a » a/Anxo and n » m. The classical 
equation of motion that we get from (23) will be exactly like the unperturbed 
one, if we replace a — > b = a + a/(47TXo) • Now, one can easily see that 
for n » m, the quantization condition 



f ° \l\V m {a)\da = {n + l/2)nh 

J 



leads to the mass quantization condition 



M(n,m) 




Vn + m + 1- <Jm + 1/2 J 



Mr 



(24) 

(25) 
0" we of 



which is slightly different from (21), (in the "limit": "m + 1/2 
course recover (21)). 

This example shows that one can get quantization condition in a more 
general case than the simple Oppenheimer-Snyder model. But if, for example, 
we make the same calculation for a minimally coupled massless scalar field, 
we get the potential 



3\ 2 



W(a) = (Snxt) 



a a — a — 



a 



k 2 



where k is a real number. But now the integral 

'\W(a)\da 



(26) 



(27) 
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diverge logarithmically, and there are no quantization conditions! The prob- 
lem is that when a — > , the potential (26) diverges. This means that the 
perturbations are not small at a — > , and we can no longer say that there is 
a collapse process. The classical solutions of (26) do not describe a collapse! 
This seems to contradict Price's results [8] that the Oppenheimer-Snyder 
model is stable. However, one should remember that to get (23) or (26), we 
neglect a very important condition, continuity at the surface of the star. So 
as a matter of fact (23) and (26) are incorrect. We can see this from the fol- 
lowing consideration: from the no-hair theorems we know that for late times 
the scalar field outside the horizon vanishes. From continuity (and the fact 
that we consider only constant scalar field inside the star), it should vanish 
also inside the star. In (26) k (and in (23) m ) is the quantum number of 
the scalar field inside the star. So if somehow (using the correct continuity 
conditions) we were able to get an effective one-dimensional Wheeler-DeWitt 
equation (like (22) with (23) (or (27))) then for late times k — > ( m — > ). 
Late times corresponds to small a, so k (and m) should be a-dependent. In 
that case (27) will not diverge. It is not easy to find the whole superspace 
with the boundary conditions [9] , and to solve the Wheeler-DeWitt equation. 
But if the Oppenheimer-Snyder model is stable, then for small perturbations, 
one should end up with an effective potential (such as we got in the over- 
simplified model of the conformally invariant scalar fielof] ) which is close to 
(14). One is tempted to make the conjecture that the WKB wave function 
of a (general) collapsing star describes a bound state, leading to quantization 
conditions. 

But one should find first the WKB solution of the full Wheeler-DeWitt equa- 
tion, or at least the full linear perturbations, which is not an easy task! 

The reason that there are no quantization condition in the standard cos- 
mological model of Friedmann universe [7] , is that in that model the potential 

6 Probably the continuity conditions will not change (24) too much, because anyway a 
and m are small and contribute a finite correction in (24). 
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is (when there is a cosmological constant^], A): 



V(a) = a 2 - H 2 a 4 

where H 2 = A/3 > 0. This potential does not describe a bound state. 
We end with some comments: 

1) It can be easily seen that our quantization conditions do not effect the 
classical collapse process, because of the correspondence limitQ 

AM M(n) - M(n - 1) 1 

« 1 . (28) 



M M(n) n 

2) We obtained a "strange" quantization condition. It depends strongly 
on the initial value Rq. As a matter of fact if R is not much greater than 
M, one should replace xl m (14) (or (24)) with 3(x /2 — sm(2xo)/4). Then 
the quantization condition (21) will be different, so the exact M is a more 
complicated function of Rq. The strange feature of our quantization con- 
dition is that the mass quanta are generally much smaller then the Planck 
mass, a more "reasonable" quantization condition should not involve Rq in 
an essential wayf]. 

3) To see if the Hawking process can be affected by the quantization 
condition, one should compare the energy of a radiated particle to the mass 
difference AM . The energy of a radiated particle is: 

1 M 2 
£? r „ ~ -ksT — ——7— . (29) 
p 2 16vrM v ; 

The Hawking process should not be affected, if E T , P » AM . In that 
case one can consider the mass as continuous. In our case this leads to 
n « Rq/Ip, which (in our geometrical units) means that 

M « Rq , (30) 



7 Otherwise the solution is not very interesting. 
8 For our sun, for example, n ~ 10 100 ! 

9 For example M ~ f(n)Mp, such that / is a function of n but not of Rq 
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consistent with reasonable cosmological initial conditions. So our quantiza- 
tion condition should not affect the Hawking process. It is interesting to 
notice that the more "reasonable" quantization condition, M(n) = f(n)Mp, 
would (very much) affect the Hawking process, because (in that case) E rp » 
AM leads to M « M P . So perhaps our "strange" quantization condition 
is not so strange? 

4) Can strong quantum effects drastically change our result? Quantum 
effects should be dominant when a — > . But if those effects should save us 
from the singularity, or make it less singular, then they must (drastically) 
increase the potential (near a = ). This means that we will have an even 
stronger binding, and quantization conditions should remain. 

Finally, it is very important to study consistently the effect of small per- 
turbations around this minisuperspace. 
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A Addition of a conformally invariant scalar 
field: 

In this Appendix we consider the addition of a conformally invariant scalar 
field, ip(t,x) . For \ < Xo ( r < r s ) we take the scalar field to be only time 
dependent, ip(t) (like in the standard minisuperspace models [7]). But for 
r > r s the scalar field must be a function of the radial coordinate too, <p(t, r) 
. The Lagrangian (8) is now: 

L = Lq + L M + = 
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ad 2 acj) 2 2 



(31) 



>r>r s 



-9 



R(l - y? 2 /6) - (Vp) s 



d 3 x 



where (pit) = a(t)ip(t) . We consider the scalar field as a small perturba- 
tion to the Oppenheimer-Snyder model. One can use "generalized Novikov 
coordinates" [10], for r > r s , whose line element is: 



ds 2 = -N 2 (R : f )df 2 + F 2 (R, i)dR 2 + r 2 (R, £)dn 2 2 
The Hamiltonian related to (31) is 

H = #l + # 2 

P, = N \—^{P 2 a -P 2 )+A'n X l{a-a Q - ( t ) 2 /a) 



(32) 



/•oo 

tf 2 = / N(R,0 



Ro 

+ { ^R 



1 



32tt 



P 2 + -P 2 + -P F P r 



F 2 Up 



(4) 



6 



Fr 2 



dF 



(33) 



where = dL/d<p , P F = dL/dF , P r = dL/df and P v = <9/<9y? . The 
corresponding Wheeler-DeWitt equation is 



where 



A 1 
A 2 



'dlT 
dN 



1 



|tf >= (A 1 +A 2 )\^ >= 



16vrxoa 



P 2 -P 2 + (87T^) 2 [a(a-G ) 



Ro 



1 



32tt 



2 P 2 + ^ + -PfA 



(34) 



(35) 



6 ) 



Fr 1 



dR 



(36) 
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and P v = —ihd/dr] , rj = (a,(j),F(R),r(R),ip(R)) . We see that A\ is a 
second order differential operator acting on a two-dimensional space. But 
A 2 acts on an infinite dimensional space (functional space). It is a very 
complicated space, with boundary conditions result from gauge fixing and 
continuity conditions at R = Rq. Because we are not going to solve the 
Wheeler-DeWitt equation in that region anyway , we will not consider all 
the details. Let z be a coordinates in that space. And let x = a , y = <f> . 
Then we make a separation of variables^ : 



V(x,y,z) = ^ 1 (x,y)q 2 (z) 



(37) 



then we get from (34) and (37) 



Ai*i(x,y) 



-a^ 2 (z) 



(38) 
(39) 



where a is a (positive) number. From (39) and (42) we get: 



B 2 r) 2 



3\ 2 



dx 2 



dy 2 



X \X — do — 



a 



*!{x,y) = 



(40) 

The solution of (40) is very simple, *$?(x, y) = ip{x)g{y). The y-part 
is just a harmonic oscillator^], for which the WKB is the exact solution. 
The frequency of this harmonic oscillator is u = IGtcXo ■ So the energy is 
huj(m + 1/2) , and from (40) we get 



ip{x) = 



(41) 



10 The coordinates x, y and z , are not independent. They are related by the continu- 
ity conditions. We will work in this simplified model and then check if our results are 
"reasonable" . 

n with a different overall sign (relative to H a )- Remember that in Wheeler's superspace, 
only the "volume coordinate" (in our case, a) is timelike, and all the other (infinite number 
of) dimensions are spacelike. 
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where 



3\ 2 



V m {x) = (Sirxti 
This is (22) and (23). 



x [ x — a — 



a 



[m+l/2)h 



(42) 
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